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Abstract 

We define the radial moment, (r), for jets produced in hadron-hadron collisions. 
It can be used as a tool for studying, as a function of the jet transverse energy and 
pseudorapidity, radiation within the jet and the quality of a perturbative description of 
the jet shape. We also discuss how non-perturbative corrections to the jet transverse 
energy affect (r). 





Prospectors for new physics congregrate at the high-energy frontier. Their searches re¬ 
quire confronting experimental data with theoretical expectations. Conhdence in their claims 
of new physics presupposes not only a proper comparison of data to theory, but also a solid 
understanding of systematic uncertainties in both theory and experiment. An interesting 
case study is the allegation of new physics based on the high-energy tail of the single-jet in¬ 
clusive transverse energy distribution measured by the CDF and D0 detectors at Fermilab 
§0. Data from both experiments agree well with both the theory (EKS for CDF and 
Jetrad iii for D0 ) and with each other for jets of transverse energy less than about 
200 GeV, while at higher energies the CDF data appear to be somewhat larger than expected 
(see for example 0]). A genuine rise above theoretical expectations at large transverse ener¬ 
gies would be the signal of one of a whole panoply of new physics possibilities. Conhdence 
in such a claim, however, requires that we hrst rule out more prosaic explanations, such 
as an uncertainty in the distributions of partons inside the colliding nucleons |§, 0, one of 
the non-perturbative inputs to the theoretical computation. In order to help resolve the 
apparent discrepancies between the experiments, and also to clarify puzzling aspects of the 
cross sections measured at lower energies |TT|, 




we feel it is important to examine other 
observable quantities in the same high transverse-energy events. For example, Ellis, Khoze 
and Stirling have proposed studying the radiation between jets in order to identify the 
underlying color structure of the hard scattering. 


In this paper, we advocate studying the radiation within a jet, and its use as a diagnostic 
tool for the quality of the perturbative description of the jet shape as a function of transverse 
energy. Both CDF and D0 have measured the transverse-energy prohle of jets, where 

the integrated density as a function of the radius r from the jet axis. 


da 


/: ET(r')dr'-^^ ^ 


/o^ ET{E)dr' 


da 

dETdrj 


In terms of particle (tower) i within the jet and lying at a distance r^jet from the jet axis, 
where, 

^liet = + (Ah^,jet)^ (1) 

we have. 


^!{r,ET,r]) = 


jets 


A'Ti(D,jet) 

The summation runs over all particles (towers) in the jet and = -^Tjet- These 

prohles have also been discussed from the theoretical point of view by Ellis, Kunzst and 
Soper DlSj IT and more recently by Klasen and Kramer . 


While one can make a qualitative or semi-quantitative comparison with theoretical pre- 
dictions,[] the language of jet prohles does not lend itself readily to the extraction of funda¬ 
mental quantities, or to identifying possible problems with underlying events or the energy 

^From the perturbative point of view suffers from large collinear logarithms for small r. 
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calibration of jets. It would be better to summarize the jet profiles in a small set of numbers, 
whose variation with Ex and rj can then be mapped more easily. The simplest such number 
is the radial moment of the transverse energy distribution within the jet, which we shall 
denote (r). 


(Si dETdri )jgtg 

(Ei (r*,jet) 


( 2 ) 


In the central region, Et ~ E, and r ^ 6] for massless partons within narrow jets, rE-r is 
then roughly 9E ~ |k| sin6', where 9 is the angle between the parton and the jet axis. This 
quantity is simply the transverse momentum of the parton with respect to the jet axis. The 
radial moment can thus be understood for narrow jets in the central region as the average 
transverse momentum with respect to the jet axis, divided by the jet’s Et- In the context 
of two-jet inclusive distributions, the triply-differential distribution (Pa/dETdriidr ]2 we have 
discussed previously extracts as much global information about the event as possible, 
from the viewpoint of perturbative QCD. The moment (r) is independent of this distribution, 
and its comparison with theoretical predictions may provide us with additional information. 
As we shall discuss below, the moment is sensitive to the size of the underlying event as well 
as to the experimental jet algorithm. 


The moment vanishes in a lowest-order calculation of jet cross sections, since jets at 
that order are approximated by lone partons; a next-to-leading order calculation of jet cross 
sections thus gives the leading non-vanishing calculation of the moment. Using the Jetrad 
implementation Q, ^ of 0{al) parton scattering processes |^, it is straightforward to 
evaluate (r) for an arbitrary jet algorithm. The dependence of the moment on Et using 
a perturbative implementation of the Snowmass algorithm with R = 0.7 is displayed in 
£g. 1. To mimic the D0 acceptance, the jets have been restricted to the central rapidity 
region, {rjl < 0.2. For reference, we use the CTEQ4M parton distribution set together with 
a renormalization/factorization scale p = 0.5Et- 


For narrow central jets, we can model the radiation pattern we expect using the collinear 
approximation to the matrix elements. We hnd, 

(r) ~ j/-^^-P(z) , (3) 

where P{z) is an Altarelli-Parisi splitting function. The dependence on the hard cross 
section, including the dependence on the parton distribution functions, disappears; only the 
dependence on remains. With 9i the angle between the jet axis and parton i, 

s = 22;(1 — 2:)£'^(1 — cos( 6 *i-I- 6 * 2 )) 

~ zil- z)E\9^ + 92Y. (4) 


Balancing of transverse momentum within the jet gives the additional relation, 

z9i = (1 - z)92. 
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Figure 1: The moment (r) as a function of the jet transverse energy at ^/s = 1800 GeV using 
the CTEQ4M parton distributions, for jets with \t]\ < 0.2 reconstructed using a perturbative 
implementation of the Snowmass algorithm (jet algorithm II dehned in eq. 7) for the jet cone 
size R = 0.7. The renormalization and factorization scales have been chosen to be equal to 
0.5Et. These parameters form our baseline choice. 


With the change of variables, ri = z9i, r 2 = {I — z)92, we have, 

(r) ~ (/ J d9i + j dz (1 - z)P{l - ^) J dd-^ . (5) 

The limits in our integral are determined by the clustering in the jet algorithm we are 
using. At next-to-leading order in perturbation theory, there are two distinct choices for 
whether or not the two partons will be clustered to form a jet of cone-size i?, depending on 
whether the parton-parton or parton-jet distance is constrained to be less than i?; 

I. ri2 < R, (6) 

II. Tjjet < R- (7) 

In this latter algorithm, two partons with equal transverse energy will combine to form a 
single jet even if ri 2 = 2R. In this case, the jet will have all the transverse energy lying on 
the edge of the jet. It is unlikely that an experimental jet, made up of transverse energy 
smeared over many calorimeter cells will reconstruct such a jet; it is much more likely to hnd 
two (smaller jets), with an overlapping region. This is impossible to model accurately at 
next-to-leading order when at most two partons can merge. However, it is usually 
approximated by the constraint, 

III. r^jet < R and ri 2 < RsepR, (8) 
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where 1 < i?sep < 2. When Rsep = 1, we find jet algorithm I, while for i?sep = 2 we obtain 
algorithm II. Even after choosing the clustering criterion, there still remains a choice of 
recombination scheme. The most common choice is the Snowmass recombination method 
^ originally designed to make semi-analytic theoretical calculations feasible. Since this 


is a reasonable approximation to the CDF and D0 recombination schemes in the central 
region^ we choose to adopt this scheme for the calculation of (r) at 0 {a^). 

For the jet cluster algorithm III, and assuming 6 i > 62 , we have the constraints 9i < R 
and 01 < (1 — z)RsepR- This translates into the bounds, 9i < R and 02 < zR/{l — z) for 
z < (i?sep —1)/-Rsep = X while 01 < (l — z)RsepR and for (i?sep —l)/.Rsep < ^ < 1/2. Summing 
over both partons (and including a factor of 2 for the case 02 > 0i), yields. 


(r)iii = dzzP{z)J^ d9i + 


TT 


+ 


dz{l — z)P{l — z) 


|•{l — z)RsepR 

dz {z)P{z) / 001 

Jo 

1 fx fzR/(l—z) '' 

002 + / dz {1 — z)P{l — z) 002 

Jo Jo 




TT 


(9) 


There are two possibilities to consider: where the parent parton is (a) a quark or (b) a gluon. 
Integrating the appropriate splitting functions, we obtain. 


V. 


Ill 

gluon 


(601n(i?3ep) + 5 + 60/i?3ep + W/Rip - S/Ri^ - 30/i?2 _ 21R,,^)N 


15 


(20 - 30 /Rsep + 30/- 20/ R^ + 6/ R^ - 3Rsep)nf 


Rluark 


(8 ln(i?sep) + 2 -|- 6/i?. 


sep 


30 

3R. 


sep 


VRipm^ - 1 ) 




( 10 ) 


Inserting the numerical values N = 3 and rif = 5 and fixing Rsep = 1 we obtain the radial 
moments for jet cluster criterion I, 


V. 


gluon 


V, 


quark 


IAN -|- Hf 
10 

3(iV2 - 1) 
4iV 


= 4.7, 

= 2 , 


( 11 ) 


while for Rsep = 2 and jet cluster criterion II, we find. 


Rfluon 


V. 


II 

quark 


(1921n(2) -43)A^ + 7n/ 
48 

( 161 n( 2 ) -3)(A^2 


= 6.36, 


8 N 


= 2.70. 


( 12 ) 


^ As detailed in ref. the D0 recombination scheme overestimates the energy of a jet in the forward 
region, and hence leads to distortions in the shape, rendering it unsuitable for the studies envisaged here. 
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Figure 2: The fraction of ‘quark’ and ‘gluon’ jets for \ri\ < 0.2 and ^/s = 1800 GeV as 
a function of the jet transverse energy for the CTEQ4M and CTEQ4HJ parton density 
functions. 


As expected, this gives rise to somewhat fatter jets. However, in both cases, Vgiuon/Vquark ~ 
2.35 so that the relative fatness of ‘gluon’ and ‘quark’ jets is preserved. Furthermore, 
~ 'Pfuark/T^luark ~ 1-35. Other choices for i?sep yield values smoothly dispersed 
between these two extremes but preserving the relative fatness of the types of jets. 


As the jet-dehning radius R increases, the collinear approximation gets worse, but none¬ 
theless we expect the general features of this simple model to survive: the moment should 
be roughly a constant times ^^(/r), where p is a scale characterizing the jet. In particular, 
it should be nearly independent of the parton distribution function. However, the parton 
distributions do play an important role. As we increase the hardness of the scattering we 
sample different types of jet. The matrix elements for the different subprocesses are (within 
color factors) rather similar for scattering at ~ 90'’, and the dominant effect is due to the 
variation of the quark and gluon parton densities with x. This is illustrated in Eg. |] where 
we show the fraction of ‘quark’ and ‘gluon’ jets at a/s = 1800 GeV as a function of the jet 
transverse energy for the GTEQ4M parton distributions [^. While ‘gluon’ jets dominate at 
Et < 120 GeV, as the quark density functions become more important, the fraction of ‘gluon’ 
jets diminishes. Different parton distributions give quite different results. In particular, the 
GTEQ fit to the supposed excess of jets high transverse momentum GTEQ4HJ containing 
an enhanced gluon at large x generates a larger fraction of ‘gluon’ jets at high Et^ 


This variation of (r) with the choice of distribution function set is contrasted in fig. ^ with 

Note that in order to fit the jet excess at high Et, the inclusive jet cross section also increased as can 
be seen in fig. ||. 
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Figure 3: (a) The ratio of (r) for the CTEQ4HJ and MRSA' {as{Mz) = 0.115) distribution 
sets to that for the CTEQ4M set, as a function of the jet transverse energy, (b) The ratio 
of the single-jet inclusive cross section computed for the CTEQ4HJ and MRSA' (as(M^) = 
0.115) distribution sets to that for the CTEQ4M set. In both cases the renormalization and 
factorization scales are equal to 0.5Et and |? 7 | < 0.2 and R = 0.7. 


the relatively greater sensitivity of the inclusive-jet spectrum. To make a fair comparison, we 
use the CTEQ4 £t with an enhanced gluon at large x that describes the single jet inclusive 
distribution well at high Ex, CTEQ4HJ and a ‘normal’ fit by the MRS collaboration 
with a similar value of as{Mz)- The relative insensitivity to parton distribution functions 
makes this moment a useful tool for measuring the strong coupling constant. 


The radial moment is also somewhat more sensitive to the jet algorithm than the single¬ 
jet inclusive cross section, as illustrated in £g. ^ for jet algorithms I and III with R = 0.7 
and Rsep = 1.3 relative to algorithm II. Once again, the jets are constrained to lie centrally, 
\t]\ < 0.2. We see that both the normalization, and to a lesser extent, the rate of decrease of 
(r) vary more than the single-jet E^ spectrum under changes of the jet algorithm. 


Of course, the calculations presented here will suffer a signihcant renormalization-scale 
dependence, as they are leading-order calculations of the radial moment; a measurement of 
as would require a next-to-leading order calculation of this quantity. To estimate the scale 
uncertainty, we compute the ratio of (r) computed with /i = 0.25Et and p = Et relative to 
the same quantity at our reference scale choice /i = 0.5Et for jet clustering algorithm II. We 
see in fig. ^ the scale variation is still sizeable and affects the normalization considerably. 
On the other hand, varying the scale hardly changes the shape. 


Armed with the perturbative predictions, we can ask whether the radial moment is sensi- 
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Figure 4: (a) The ratio of (r) for the jet algorithms defined in the text as function of 
the jet transverse energy, (b) The ratio of the single-jet inclusive cross section for the 
same jet algorithms. In both cases, CTEQ4M parton distributions have been used with 
renormalization and factorization scales equal to 0.5Er and |? 7 | < 0.2 and R = 0.7. 


tive to non-perturbative contributions, either to the jet energy or to the pattern of radiation 
within the jet. These may arise from various sources, such as power (1/Q) corrections to 
the hard scattering process, spectator interactions, or soft hadrons in the proton/antiproton 
remnants spilling into the jet cone. Another source of uncertainty is detector-related effects 
such as the jet energy-scale uncertainty. For example, the debris of the hard scattering tends 
to form a roughly uniform distribution in pseudorapidity and azimuth. Shower models are 
used to estimate this underlying event and thereby correct the jet energy. These have been 
extensively tuned at moderate energies where the data are plentiful. However, if the under¬ 
lying event was not properly modelled as a function of Et, a residual jet energy correction 
could occur. 


In the simplest string-like hadronization model p4 
order C 


we expect to hnd a contribution of 


1 GeV to the transverse momentum pt with respect to the jet axis in a central 
jet. Since the average pt is roughly {r)ET, this implies a non-perturbative correction of order 
C/{{r)E'j')- We expect the corrections to the radial moment to be of the same order. The 
underlying event is also expected to contribute a. 1/Q correction. If we assume that the 
energy is deposited uniformly, we expect a contribution to (r) of order 


t^Eu 

Ej' 



(13) 


where Eu is the underlying event transverse energy density (per unit rapidity per radian). 
In minimum-bias events, this is ~ 0.55 ±0.1 GeV according to ref. m- Some of this is 
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Figure 5: The radial moment (r) for renormalization/factorization scale fi = XEt relative 
to that for fi = Et for central jets, \r]\ < 0.2 using the cluster criterion II and the Snowmass 
recombination procedure. The solid (dashed) lines shown the change for A = 0.5 (2.0). In 
all cases, the MRSA' parton distributions have been utilized. 


subtracted as part of the experimental analysis. Since (r) varies from ~ 0.25 to ~ 0.07, this 
correction should anyway be much smaller than that due to hadronization. 

A determination of «*, which is less sensitive to the parton distributions than other 
determinations at hadron colliders, and the possible identihcation of detector and/or physics 
based non-perturbative contributions are both interesting to pursue. This gives the jet-shape 
analysis a well dehned physics goal, and motivates the determination of jet shapes over a 
large range of transverse energies. While a sensible analysis requires at least the next-to- 
leading order contributions to the jet shape, we can already confront the published CDF 
and D0 data on the quantity T(r, Et, rf) dehned in eq. 1 with the leading-order calculations. 
This allows us to explore the potential utility of this measurement, and indicate expected 
uncertainties and problems. 

The current published data on the jet shape from both CDF and D0 are very limited. 
The CDF collaboration published results based on 4.2 pb“^ with only three bins of transverse 
energy, while D0 used 13 pb“^ with four bins in transverse energy. A full analysis of current 
data would include on the order of 100 pb“^ from each collaboration, preferably using the 
same binning in Et as used in the one-jet inclusive transverse-energy distribution. Substan¬ 
tial experimental improvements to the results presented here are thus possible. In addition, 
we expect the next-to-leading order theoretical corrections to be available soon. The lat¬ 
ter will enable us to extract the next-to-leading order as and determine the uncertainty on 
the radial moment, and will be crucial in understanding the magnitude of non-perturbative 





^T-bin (GeV) 

{Et) (GeV) 

experiment 

experiment 

40-60 

45 

GDF 

0.31 ± 0.02 

45-70 

53 

D0 

0.29 ± 0.01 

65-90 

75 

GDF 

0.24 ± 0.02 

70-105 

81 

D0 

0.23 ± 0.01 

95-120 

100 

GDF 

0.22 ± 0.02 

105-140 

118 

D0 

0.19 ± 0.01 

140-900 

166 

D0 

0.17 ± 0.01 


Table 1: The D0 and CDF experimental resnlts for (r) extracted from the pnblished jet 
shapes. 


effects. 


We mnst hrst determine the radial moment (r) from the pnblished resnlts for the qnantity 
T. In an analysis along the lines snggested in this paper one wonld of course determine the 
radial moment as a function of transverse energy directly from the data without resorting 
to extracting the explicit shapes. This would reduce the experimental uncertainties signif¬ 
icantly, especially the systematic uncertainties. For now, we have to reconstruct the radial 
moment from the quantity T; 


{r{ET,r])) 



d^{r,ET,r]) 

r --- dr 

dr 


(14) 


The values of (r) extracted from the published data are shown in Table 1. Both experiments 
use R = 1, but the allowed range of rapidity for the jets is different. CDF consider jets in 
the range 0.1 < \ri\ < 0.7 while D0 have a tighter restiction, \ri\ < 0.2. 


For the theoretical prediction, we could repeat the experimental extraction, by hrst de¬ 
termining T using the same binning as the experiments and extract (r) from this quantity. 
This will simulate any effects on (r) due to binning and functional parametrization. How¬ 
ever, these effects are in practice small and instead we determine (r) directly as was done 
in £g. ID, For the theoretical prediction, we use the CTEQ4M parton distribution set, with 
as{Mz) = 0.116. The renormalization and factorization scales were chosen to be half the 
jet transverse energy. To match the experimental cuts, we use R = 1 and i?sep = 1-3 was 
used to simulate the jet splitting and merging effects, which are not modelled at this order 
in perturbation theory. 

The extracted experimental values for (r) are shown in £g. ^ along with the theoretical 
prediction as a function of Et- We see that the CDF and D0 points are compatible with 
each other (despite the fact that the experiments use slightly different rapidity cuts on the 
jet). On the other hand, the leading-order predictions cannot describe the data, as its Et 
dependence is substantially different from that of the data. As was already demonstrated 


9 










Figure 6: The extracted values for (r) for D0 (open circles) and CDF (solid circles). The 
theoretical prediction for |p| < 0.2, CT 4M parton distributions, R = 1 and i?sep = 1-3 
and p = 0.5Et is shown as a solid line. The fitted theoretical result with an additional 
non-perturbative component as described in eq. and C = 1.09 and K = 0.56 is given by 
the dashed line; the dotted lines indicate the la variations in the £t. 


in fig. 1^ while the normalization of the leading-order prediction has a strong dependence 
on the choice of the renormalization/factorization scale, the shape is mostly independent 
of this choice. Thus, although we assume that higher-order corrections will modify the 
normalization, it is reasonable to assume that they will not modify the shape substantially. 

This suggests that the explanation of the shape of the data lies in hadronization effects 
(i.e. power corrections). To test this hypothesis, we will assume that the perturbative 
prediction (r)perturbative is a constant K times the leading-order prediction {r)Lo- We model 
the hadronization effects as described earlier, so that 

T'') hadronic 


where the non-perturbative scale C should be close to 1 GeV. We then perform a fit of 
this model to the data varying both the iF-factor and the hadronization scale C. The fit 
result with a hadronization scale C = 1.09 ± 0.05 GeV and K = 0.56 ± 0.11 is shown in 
fig. 1^. Other choices of renormalization scale give roughly the same value of C. This is not 
surprising because C changes the shape of the curve while the renormalization scale affects 
primarily the normalization. This shows that the above model might explain the observed 


ij ') partonic T 

^(r)LO ' 


(J 


lr-\ pF* 

y /partomc^T 


c 




(15) 
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differences between the leading-order predictions and the data. The it'-factor has a strong 
dependence on the renormalization scale, as expected. 


Only a next-to-leading order calculation of (r) can conhrm that power corrections, rather 
than higher-order corrections, are responsible for the observed behavior of the data. Such a 
calculation requires the one-loop hve-parton matrix elements as well as the six-parton 
tree-level ones all of which are available in the literature. With such NLO predictions, 
one could extract both an NLO as (sensitive to the normalization and the high-T^'r depen¬ 
dence of (r)) and the hadronization scale C (sensitive to the Iow-Et dependence of (r)). In 
order to perform such an extraction, however, we need a much more extensive measurement 
of the Et dependence of (r). With current CDF and D0 data sets, it should be possible to 
cover a much larger E^r^^ range (both lower and higher than the current published results) 
with a hner binning. This would enable a study of power corrections and their impact on 
the uncertainty in an extracted as- 

A few remarks are in order. The current iterative cone algorithm will be rather problem¬ 
atic in such an NLO calculation (these are the problems one would encounter in a next-to- 
next-to-leading order calculation of jet differential cross sections). For the studies performed 
in this paper, we used a phenomenological parameter i?gep = 1-3 to model jet splitting and 
merging. At the next order in perturbation theory, one should either abandon this param¬ 
eter (relying on the theory to model jet splitting and merging), or else one must introduce 
an additional For the purposes considered here, i?sep is anyway an especially suspect 

parameter because it is purely phenomenological and may well be Et dependent [ 0 - In¬ 
deed, as we have seen in this paper, jets get narrower as their energy increases. With this 
observation one could argue that i?sep should in fact decrease with increasing energy of the 
jet. A simple way out of this morass is to use the Kt jet algorithm pR [1^. It has good 


perturbative behaviour and no additional phenomenological parameters are needed. 


The analysis discussed in the present paper suggests that it should be possible to use the 
hadronic dijet system to determine both as and the proton’s gluon distribution with great 
accuracy. The gluon distribution would be determined by the triply-differential distribution, 
and, once the next-to-leading order calculations are available, a* would be determined by the 
radial moment. Only the quark distributions would be taken from deeply-inelastic scattering 
data. Use of the full current data set would allow a much hner jet Et binning as well as much 
higher jet Et to be used, reducing the experimental uncertainties on (r) signihcantly reduced, 
and thereby reducing the signihcance of power corrections. If the higher-order perturbative 
corrections are also small, as suggested by the phenomenological study performed in this 
paper, the theoretical uncertainties should be competitively small as well. We may also 
expect the high-U^ jet prohles to be especially sensitive to jet production from new physics, 
and thus the Et dependence of the radial moment should be a good probe of the presence 
(or absence) of physics beyond the standard model. 
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